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Abstract 

Binding energy of symmetric nuclear matter can be accessed straightfor- 
wardly with the textbook mass-formula and a sample of nuclear masses. We 
show that, with a minimally modified formula (along the lines of the droplet 
model), the symmetry energy of nuclear matter can be accessed nearly as 
easily. Elementary considerations for a macroscopic nucleus show that the 
surface tension needs to depend on asymmetry. That dependence modifies 
the surface energy and implies the emergence of asymmetry skin. In the mass 
formula, the volume and surface and (a)symmetry energies combine as ener- 
gies of two connected capacitors, with the volume and surface capacitances 
proportional to the volume and area, respectively. The net asymmetry parti- 
tions itself into volume and surface contributions in proportion to the capaci- 
tances. A combination of data on skin sizes and masses constrains the volume 
symmetry parameter to 27 MeV < a < 31 MeV and the volume-to-surface 
symmetry-parameter ratio to 2.0 < a/ (5 < 2.8. In Thomas- Fermi theory, the 
surface asymmetry-capacitance stems from a drop of the symmetry energy 
per nucleon S with density. We establish limits on the drop at half of normal 
density, to 0.57 < S(po/2)/ S(po) < 0.83. In considering the feeding of sur- 
face by an asymmetry flux from interior, we obtain a universal condition for 
the collective asymmetry oscillations, in terms of the asymmetry-capacitance 
ratio. 
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I. INTRODUCTION 



Changes in nuclear properties and in reaction dynamics under the changing nucleon con- 
tent of a system attract attention in the context of exotic nuclear beams [1]. Determination 
of the density dependence of the symmetry energy, describing rise in the nuclear energy with 
an increasing asymmetry between neutrons and protons, could permit extrapolations to the 
limit of neutron matter and determinations of neutron-star properties [1-3]. Central nuclear 
reactions can provide varying nuclear densities [1] but constraints on the symmetry energy 
from structure, pertaining to normal and subnormal densities, are also very important. As 
an example, the nuclei with a significant excess of neutrons over protons tend to be charac- 
terized by larger radii for the neutron than proton distributions, i.e. exhibit neutron skins. 
It has been noted, within nuclear mean-field calculations [4-7] , that the skin size could yield 
information on the derivative of symmetry energy with respect to density. For knowledge 
to be easily transferrable between the structure and the reactions, which bring in their own 
complexity, it is important to understand the role of symmetry energy in the structure at a 
basic level. For a utility in the reactions, importantly, any quantitative conclusions from the 
structure cannot be just tied to a specific structure model but need to pertain to universal 
properties of the medium, even if coarsely. 

Most elementary information on nuclear properties stems from the semiempirical energy 
formula describing the average changes in nuclear binding energy with the nucleon content 
[8,9]. The volume and surface terms in the standard formula pertain to symmetric systems. 
The volume coefficient provides the binding energy per nucleon in the symmetric nuclear 
matter; the surface coefficient, up to a factor, provides the surface tension. Notably, the 
symmetry term in the standard energy formula has a volume character only. However, as 
we point out in this paper, basic intrinsic consistency of the semiempirical formula requires 
the inclusion of a surface symmetry term. Quantitatively, such a term turns out further 
to be essential for describing energies of light asymmetric nuclei. Moreover, inclusion of 
such a term is essential for drawing conclusions about the symmetry energy in general, from 
the nuclear structure. The surface symmetry energy appears in the literature in advanced 
energy formulas, generally in two different type of forms, in an expanded form such as in the 
liquid drop model by Myers and Swiatecki [10-13] or in more complex forms of the liquid 
droplet model by the above authors [14-16,12]. Both formulations are currently employed 
in the descriptions of nuclei [12,13,17]. 

In this work we follow a general type of considerations for the energy of a nucleus treated 
as a binary macroscopic system with a surface. Similar type of considerations can apply 
to other binary systems with an interface, such as between two phases [18]. We show that, 
when the surface energy is affected by the particle asymmetry within the system, then the 
thermodynamic consistency requires that some of the asymmetry moves to the surface, i.e. 
an asymmetry skin develops. Minimization of the net nuclear energy with respect to the 
partition of asymmetry produces an expression for the symmetry energy such as in the more 
recent versions of the droplet model [15,16]. (In the droplet model, the skin size is a basic 
parameter and, as such, one of the starting points.) We show that the potentially confusing 
expression for the energy is easily comprehended using a capacitor analogy. The symmetry 
energy, quadratic in the asymmetry, can be interpreted in terms of the energy of a capacitor 
for asymmetry, with the capacitance for asymmetry proportional to the nuclear volume in 
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the standard energy formula [8,9]. However, in the actual nuclei, the surface provides an 
additional capacitance for the asymmetry, proportional to the surface area. In the end, the 
net symmetry energy is given in terms of the net asymmetry squared divided by the net 
capacitance, where one term in the capacitance grows as the volume and the other as the 
area. Notably, the liquid-drop type formula follows when the surface capacitance is small 
compared to the volume capacitance. In general, when ignoring the surface term in the 
nuclear energy formula, the whole asymmetry capacitance is attributed to the interior and, 
thus, the volume capacitance is overestimated which is equivalent to underestimating the 
volume symmetry coefficient. 

On adding a single parameter to the textbook nuclear energy formula, representing the 
surface symmetry energy, we fit the measured ground-state energies [19] to establish limits 
on the symmetry energy coefficients. Notably, the tradition within the description of nu- 
clear energies in terms of advanced formulas [10-16,12], often supplemented with structure 
calculations, has been to quote optimal or preferred values for utilized constants. While 
this is fine as far as parameters of a model are concerned, this is not sufficient when the 
fundamental matter properties are being sought. As a drastic example, the authors of [12] 
produce parameters in two different descriptions of nuclear energies that yield quite different 
surface symmetry properties for nuclear matter, seemingly without a concern. For the fun- 
damental properties, the validity of employed theory needs to be assessed and, importantly, 
the features of the optimization, i.e. the value for the property needs to be quoted together 
with an uncertainty. However, this is normally not done in the context of advanced mass 
formulas. Obviously, if the fit minimum is very flat, or nearly degenerate, the property is 
not well constrained. One plain example of how things can go wrong with a fit is when a 
parameter is included that has little bearing on the quality of the fit, yet the optimization 
produces a value for the parameter. For validity of conclusions, the importance of the spe- 
cific parameter for the fit needs to be demonstrated. For a multitude of parameters, the 
issues of fit stability may need to be explored. 

In analyzing our fit, we demonstrate that the surface symmetry term is required, in 
order to describe the energies of light asymmetric nuclei at the level similar to the other 
nuclei. However, once the surface symmetry term is included, the minimum for the fit to 
masses is relatively flat when the surface and volume symmetry coefficients are varied along 
a correlation valley enforced by the need to reproduce, for heavy nuclei, the magnitude of 
symmetry coefficient in the standard formula. 

We subsequently turn to the analysis of the size of asymmetry skins. The dependence 
of the size on net nuclear asymmetry and on mass probes the ratio of volume to surface 
symmetry coefficients. This is because the size represents the surface asymmetry charge 
which, in relation to the net asymmetry charge, tests the ratio of surface to volume capaci- 
tances. The dependence of skin sizes on the difference of proton-neutron separation energies 
probes, on the other hand, the surface energy coefficient. This is because, in the capacitor 
analogy, the difference of separation energies represents the voltage and the dependence of 
surface asymmetry charge on voltage tests the surface capacitance. On analyzing a variety 
of data on skin sizes [20-27], upon inclusion of Coulomb corrections, we obtain constraints 
that, in combination with the constraints from mass analysis, produce a relatively narrow 
correlation range for the volume and surface symmetry coefficients. 

In illustrating further the importance of the surface symmetry energy for nuclei, we de- 
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rive new solutions for asymmetry oscillations in nuclei, with the asymmetry flux feeding the 
surface. This extends the treatment of isoscalar modes where the internal compression and 
the surface deformation are coupled [28]. We obtain a boundary condition of a simple univer- 
sal form, selecting wavenumbers for different multipolarities of vibrations. We find that the 
use of a realistic surface symmetry coefficient improves significantly the description of the 
dependence of giant dipole resonance energy on mass, as compared to the use of assumption 
that the surface accepts no asymmetry. The mode of oscillation, in our description, changes 
from one that is closer to the Goldhaber- Teller (i.e. surface) type [29] to one closer to the 
Steinwedel- Jensen (i.e. volume) type [30] as the nuclear mass increases. 

As a final step, we analyze the nature of the surface symmetry energy at a microscopic 
level within the extended Thomas-Fermi model. We show that an extra capacitance for the 
asymmetry associated with the surface is due to the drop of symmetry energy per nucleon 
with density (lower energy penalty = higher capacitance). If the symmetry energy per 
nucleon is density independent, no extra surface capacitance emerges. Within the Thomas- 
Fermi approximation, for a given density profile of symmetric matter, the ratio of the volume 
to surface energy parameters is solely related to the shape of the density dependence of 
symmetry energy. Equipped with this knowledge, we turn to an analysis of the mean-field 
structure calculations in the literature [5-7], where numerical experimentation was used to 
find a relation between the skin sizes and features of the nuclear equation of state. The 
combination of our earlier results and those in the Thomas-Fermi approximation indicates 
that the skin for a given nucleus is primarily related to the shape of the density dependence 
of symmetry energy per nucleon. Using the results from the structure calculations [7], we are 
able to narrow the bracket on the density dependence of symmetry energy by nearly a factor 
of 2 compared to what was done before [7]. Also, within the Thomas- Fermi approximation, 
we cross-test the validity of our general results for skins. 

The paper is organized as follows. In the next section, we discuss the symmetry energy 
within the binding-energy formula and we analyze known masses in terms of the modified 
formula, to constrain the symmetry energy parameters. In Sec. Ill, we derive expressions 
for the nuclear skin size in terms of the symmetry energy parameters, as a function of 
mass and charge numbers and as a function of separation energies. We analyze data in 
terms of obtained expressions to get additional constraints on the symmetry parameters, 
narrowing the range allowed by data. The asymmetry oscillations are discussed in Sec. IV. 
The constraints on the low-density behavior of symmetry energy are discussed in Sec. V. 
Our results are recapitulated in Sec. VI. 



We start out with a reminder of elementary issues. The Weizsacker nuclear energy 
formula [8,9] 



II. BINDING ENERGY FORMULA 



A. Elementary Considerations 



E = -a v A + a s A 2/3 + a, 



a 



{N-Zf 
A 



+ ac 



Z(Z-l) 
A 1 / 3 



(1) 
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separates out different contributions to the nuclear energy. The r.h.s. negative volume- 
term represents the energy in the bulk limit of symmetric nuclear-matter in the absence 
of Coulomb interactions. Proportionality to the mass number A for this largest in magni- 
tude term, together with an approximately constant nucleon density of p = 3/(47rrQ) ~ 
0.160 fm~ 3 in the nuclear interiors, underscore the short-range nature of nuclear interactions. 
The three subsequent terms in the Weizsacker formula represent an increase in the energy 
compared to the bulk symmetric limit. These terms are associated, respectively, with the 
size of the surface, with an asymmetry in the number of protons Z and neutrons N and 
with the Coulomb repulsion. Regarding the (a) symmetry term in the Weizsacker formula, 
we observe that this term scales like A when the nucleon number is changed and that we 
do not have a surface (a)symmetry energy term, scaling like A 2 / 3 . The N-Z asymmetry 
reduces the binding because it enhances the role of Pauli principle and because the like- 
nucleon interactions are less attractive than the neutron-proton interactions. The precise 
symmetry of that term around N — Z = expresses the symmetry of nuclear forces under 
the neutron-proton exchange, or charge symmetry. 

With the radius ro of the nuclear volume per nucleon, the size of the nuclear surface is 
S = An Tq A 2 / 3 , and we can write the surface energy as 

E s = a s A 2 ' 3 = -^- 2 S. (2) 
Air Tq 

The tension a is the derivative of energy with respect to surface area: 

dE a s 



dS Anr 2 



(3) 



The tension represents a work per unit area needed to create the surface. The work is 
required because the nucleons at the surface are less bound than in the interior. However, 
we know that the N-Z asymmetry reduces the interior binding. Thus, with an increasing 
asymmetry, the work required per unit surface, the tension, should drop. Coulomb effects 
will also affect the surface, but will be dealt with later. 

As a microscopic quantity, the tension should depend on a microscopic quantity associ- 
ated with the particle asymmetry. The respective quantity common for different subsystems 
in contact is the chemical potential. With the nucleon chemical potentials given by 

fin= dN and fi > = dZ' () 

the potential conjugate to the asymmetry N—Z is half of the the nucleon potential difference, 
fJ>a = (fJ> n ~ / i p)/2, and the change in energy associated with the particle numbers is: 

dE = \i n dN + \i p dZ = jj, dA + /i a d(N — Z) , (5) 

with the median chemical potential ji = (fx n + /i p )/2. 

Given that the tension should drop no matter, under the charge symmetry, whether a 
neutron or proton excess develops, the tension needs to be, in the lowest order, quadratic in 
the asymmetry potential \i a : 

a = a --ffi 2 a (6) 
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where 7 > 0. As the tension is a derivative of the surface energy, the dependence of a on 
fi a implies the dependence of the surface energy on fi a , with serious consequences. 

To demonstrate these consequences, let us construct the thermodynamic potential gen- 
erated by a Legendre transformation from the asymmetry N — Z to the potential /i a : 

$ = f ji a (N-Z)-E. (7) 

For definiteness, we next consider changes in the net system when the asymmetry changes, 
while A and shape remain fixed. The derivative of <3> with respect to /i a gives the net 
asymmetry for the nucleus N — Z: 

,r r, d(N-Z) BE Ar „ 
j— = N-Z + fi a ^— L- =N-Z, (8) 

cf. (5). On the r.h.s. of (7), we can separate the energy E into the surface Eg and volume Ey 
contributions, E = E$ + Ey. When the asymmetry changes, the change in volume energy 
is dEy = n a d(Ny — Zy), where d(Ny — Zy) is the change in particle asymmetry within the 
nuclear volume, expected equal to d(N — Z). If we now substitute (7) into the l.h.s. of (8) 
and compare to the r.h.s. there, we find 

V 0/J, a d/jL a J d/J, a 

Here we can see, though, that we will run into a contradiction if we identify changes in the 
net nuclear asymmetry N — Z with the changes in the interior asymmetry Ny — Zy, when 
the surface energy E$ depends on chemical potential fi a . The net and interior asymmetries 
must be different, and the difference N$ — Z$ = N — Z — (Ny — Zy) must, as a consequence, 
reside within the nuclear surface. To tackle the conceptual issue of how a surface can contain 
a particle number, we must consider how the surface quantities get to be defined. 

For an interface between phases in equilibrium, Gibbs [18] suggested to take a small but 
macroscopic volume surrounding an element of the interface and, besides the actual system, 
consider an idealized reference system, cf. Fig. 1. In the idealized system, the densities of 
thermodynamic quantities, such as the particle density, change abruptly at the interface 
from a value characterizing one of the phases away from the interface to the value for the 
other phase. The surface share F s for the surface element S is then the difference between 
the value of the quantity for an actual system F and the value for an idealized system F^: 
Fs = F — Fid- In terms of densities, for the case of a nucleus in equilibrium with the vacuum, 
this further yields for a quantity F: Fs = f S = F — V fy, where fy is the density of F 
characteristic for the interior and f represents the surface density. The problem is, however, 
that the result for the surface depends on a precise position assumed for the interface, and 
a condition may be imposed further to specify the position. For a nucleus it is natural 
to assume that the surface contains no mass number and, for a spherical nucleus, this is 
equivalent to taking the equivalent sharp-sphere radius for the surface of R = r A 1 / 3 . 

Only one condition related to the particle numbers, however, suffices to localize the 
surface. In a two-component system, the attempts to localize the surface using one of the 
components or the other would lead to two different positions. If the surface is localized 
using the net particle density, a finite surface asymmetry Ns — Zs results, even though the 



6 



surface has a vanishing net nucleon number As = N$ + Z$ = 0. For a spherical nucleus, 
this corresponds to having two different equivalent sharp-sphere radii for neutrons R n and 
protons R p , bracketing the sharp-sphere radius for the matter R, cf. Fig. 2. 

The Gibbs prescription can be used for clarifications, but also as a direct tool, as e.g. later 
illustrated by a consideration in Subsection V A or when pursuing curvature effects. We next 
relate the surface asymmetry excess to the asymmetry potential. Either by following the 
Gibbs prescription, with a small volume surrounding some surface element, or by isolating 
contributions pertaining to the surface element in the energy considerations before, we find 

dE s = adS + n a d(N s -Z s ), (10) 

under independent variations of asymmetry and of the surface element size. From changes 
in Es, when taking an increasing fraction of the net surface, we find that 

E s = aS + pi a (Ns-Zs). (11) 

On taking a differential of the above, with the result 

dE s = S da + a dS + (N s - Z s ) dfi a + ^ d(N s - Z s ) , (12) 

and on subtracting from this result the relation (10), side by side, we obtain a Gibbs-Duhem 
type relation. That last relation yields 

Ns-Zs da 

2 7 A*a, (13) 



dji a 



where, at the end, we employed Eq. (6). 
With (13) and (6), we get from (11): 



l0 , 1 (N s -z s y 



E s = a S + lf , 2 a S = E s + - " - '- , (14) 

where E s is the surface energy in absence of the surface asymmetry. At this point, we may 
cross-check the consistency of our result (14) with its derivation. Thus, from (14), we find 



dE s = dE s + - - ^ » d{N s - Z S ) - »' dS 

= (a - 7 nl) dS + n a d(N s -Z s ) = adS + » a d(N s - Z s ) , (15) 



which agrees with (10) and where we utilized (13) and dE s = a dS. In connection with 
(14), we may note that although the surface tension decreases with the asymmetry, Eq. (6), 
the net surface energy increases with the asymmetry. 

On account of charge symmetry for the interior nuclear energy per unit volume, to lowest 
order in the asymmetry, the volume energy will have the form 

E V = E° V + Vap (^A 2 = K + a {Ny - Zv) \ (16) 



where we factorized the coefficient in the energy density so as to get the net volume energy 
in conformance with (1) and where Ey (= ay A) is the energy for a symmetric system. 
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For an analogy with (16), we introduce a constant with the dimension of energy in Eg, 
{3 = 1/(16 7T7q 7), to write the r.h.s. of (14) as 

Es = ^P^£ l - (IT) 

In the volume energy we use a constant a in place of a a of (1) to make a distinction between 
the symmetry constant in the bulk and the constant in use in the binding-energy formula. 
As will be next seen, a a in use generally results from a combination of a and (3. 
Minimization of the net system energy E under the constraint of net asymmetry, 

E = E V + E S with N-Z = N V -Z V + N S -Z S , (18) 

yields the condition of equal chemical asymmetry potentials for the interior and surface, 

dE _ dE v dE s 

d(N v - Z v ) ~ d{N v - Z v ) ~~ d(N s - Z s ) 

. 2a ^_ 2 ^. (19) 

With this, we find for the asymmetry ratios 

N S ~ Z S _ a N S - Z S _ Ny-Zy _ 1 / 9f| x 

N v -Z v /3AW N-Z l + #^ 1 + ^AW 1 } 

The net energy then becomes equal to 

a (N-Z) 2 (3 (N-Z) 2 







E — Ey + Ea — En + — — — ■ - - -7— + 



A (1 + I A"V3)2 ^2/3 + £ A1/3 y 

F 1 g (Nj-Zf 



B. Capacitors for Asymmetry 

If a problem can be mapped onto one that is well understood, existing intuitions can be 
exploited to make predictions and lay out results with little need for any derivations, and 
this pertains to the last results. The quadratic energies (16) and (17), with the constraint 
(18), invite e.g. an analogy to two connected springs, of spring constants ^ and , 
respectively, stretched between separated supports. The context, however, suggests an even 
more appropriate analogy to two connected capacitors carrying charges qx = Nx — Zx 
(X = V, S), of capacitances Cy = A. ari d Cs = 4^-, respectively, to yield capacitor energies 

2 

Ex = E x + ^fp. As the asymmetry moves in-between the capacitors, the equilibrium 
situation is that of a charge in proportion to capacitance, yielding 

N s - Z s _ C s a N s -Z s C s I , . 

Ny-Zy Cy (3 A 1 / 3 N-Z Cy + C S 1 + - A 1 / 3 ' 1 ' 



8 



The net capacitor energy, finally, is given in terms of the net charge and net capacitance 



p_ F o | g 2 _ F o | (N-ZY _ a (N- Zf 



a 



The connected capacitor and the nuclear symmetry energy problems are, at this level, com- 
pletely equivalent. 

The symmetry energy of (23) should, generally, replace that in the energy formula (1), 
with (3 — > oo (7 — > 0, in (6)) obviously corresponding to the standard formula. In (22), 
we observe that the surface fraction of the asymmetry excess decreases as nuclear mass 
increases, because the surface capacitance grows only as A 2 / 3 . It may be tempting to expand 
the symmetry energy (23) in the surface-to-volume capacitance fraction. However, whether 
this can be done depends on the volume-to-surface coefficient ratio a//3. For large enough 
ratios and small A, we, in fact, principally might be able to ignore 1 in the denominator 
and get the symmetry energy such as characteristic for the surface storage of asymmetry. 
(This would express the fact that, just as the volume capacitance grows faster than surface 
capacitance, it also drops faster.) Later in this and in the following section we shall see what 
a /(5 values are favored by data. 



C. Coulomb Effects 

So far, we have not considered the effect of Coulomb interactions on the surface asym- 
metry and, generally, the interplay of Coulomb and symmetry energies in the net energy. 
The Coulomb energy may be represented as 



e 2 1 Ci „o „ „ 1 

"'-^ 3 r v.-, 



En = E?r + Ep/o + Eq = — ( - Z?/ + Zv Zq + — Zl 



a C ( r?2 , 5 ry ry ,5 



A i/3 y Z v + ^Z V Z S + -Z 2 j (24) 

where the three terms correspond to the volume, volume-surface and surface interactions, 
respectively, and Zy = Z — Zg. 

The combination of the symmetry energies (16) and (17) and the Coulomb energy (24) 
is quadratic in the surface excess and is easily found to minimize at: 

N-Z-^ZA 2 / 3 N-Z-^ZA 2 ' 3 

1+g v:v^ a ■ (25) 

a 48 a a 

The last approximation expresses the fact the Coulomb interactions represent a very soft 
spring in the energy combination, compared to the volume asymmetry spring. Compared 
to (22), we see that, as might be expected, the Coulomb interactions try to tilt the surface 
asymmetry towards a proton excess. Quantitatively, however, the Coulomb correction to 
the asymmetry excess is quite negligible for light and medium nuclei and acquires only a 
limited significance for heavy nuclei. 

As to the energies, even when the Coulomb correction to the surface asymmetry has a 
significance, the corresponding changes in the asymmetry energy are quite ignorable. This 
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because around the minimum, specified by (22), the changes in the asymmetry energy are 
quadratic in the asymmetry correction. On the other hand, the correction to the Coulomb 
energy, from the development of surface asymmetry, 



AE C ^ 



dZ, 



a c Z 2 N-Z 1 

Zr ^lTi + ^i/3' (26) 



may have a moderate significance within the net energy and, in heavy nuclei, can reach a 
comparable magnitude to the correction to the Coulomb energy from surface diffuseness. 



D. Relation to the Droplet Model 

The droplet model [14] utilizes the relative displacement of neutron and proton surfaces 
as one of the different parameters describing a nucleus. A global energy functional is con- 
structed for a nucleus, incorporating an effect of the surface displacement on the net energy. 
This functional is minimized with respect to the different parameters in order to arrive at 
the ground-state energy. 

While the droplet model contains a variety of potentially useful details, these can easily 
shroud the issues of symmetry energy. In fact, our own understanding of that aspect of the 
model stems largely from a mapping of the results there onto those in the present paper. 
In particular, if the final result for energy in the droplet model in [15] (and in later papers 
on the model) is expanded in the net particle asymmetry, then a symmetry energy emerges 
such as in Eq. (21), provided our constants are identified with those in the droplet model 
as: 



4 §H 
9 Q ~ T^TP/J 



a = J and f3 = = - — . (27) 



The individual constants H, P and G, found equal to G — 3 J P/2Q, describe the assumed 
dependence [14] of the surface energy on bulk asymmetry and on normalized size of the asym- 
metry skin. The different definitions of the surface energy within the model are discussed 
in Refs. [31,32]. 

If the net energy in the droplet model [15,12] is expanded simultaneously in terms of the 
net asymmetry and in the Coulomb coefficient ac, then a correction to the energy is found 
in the form such as in (26). 



E. Comparison to Binding-Energy Data 

We next examine what can be learnt on the nuclear symmetry energy from measured 
binding energies [19], B = —E, employing the results from above, by following an energy 
formula with only one more adjustable parameter than in the textbooks [28], to represent 
the surface symmetry energy. Specifically, we fit the measured energies with 

,2/3 » (N-Zf Z(Z-l) 

E=-a r A + a s A^ + TTfJ ^^ — +°c Al/ l il + l } -S. (28) 



10 



Compared to (1), we changed the symmetry term to (23) and we added the standard pairing 
term 5 = ia^A" 1 / 2 ^. We either carry out fits with A = in the Coulomb term, i.e. the 
energy formula modified, from the standard, only by the inclusion of surface symmetry 
energy, or with 

57T 2 _jf 1 N-Z 

6 rlA^ ~ 1 + &A 1 /* QZ ' ( ' 

including the effects of surface diffuseness and of the asymmetry skin (26) on the Coulomb 
energy. In the first r.h.s. term above, d ~ 0.55 fm is the diffuseness parameter in the Fermi 
function from the parametrization of nuclear charge distributions [33], 

P(r) = l + exp[(r-c)/rf] ; (30) 

that term stems from an expansion of the Coulomb energy for (30) in terms of d. 

Qualitatively, results from the two fits indicated above are similar but they differ in 
quantitative details. Obviously, inclusion of the volume and surface competition, in the 
symmetry energy, changes the dependence of the net energy on asymmetry as a function 
of mass. However, the Coulomb energy is also a source of the dependence on asymmetry 
and deciding on the details in one energy contribution depends on the details in the other. 
We believe that more reliable results are obtained with the finite A, Eq. (29), in (28) as we 
only incorporate information available from electron scattering or following from intrinsic 
consistency of the considerations. The number of independent parameters in the fit with 
(29) is the same as in the A = fit, precluding, in practice, a fit runaway in the parameter 
space. 

In the fits, we optimize the sum of absolute deviations of the theoretical energies from 
experimental, rather than the sum of the deviations squared, as the intention is to provide 
an average description of the nuclear features. Results, however, weakly depend on that 
choice. In the optimization, we include the energies of all nuclei with A > 2, but exclude 
the energies deduced from trends [19]. We either let all parameters in (28) vary in an 
unconstrained manner or we constrain a/ (3 or both of the symmetry parameters. 

The average per nucleus deviation of the theoretical energy from experimental exhibits 
a narrow valley in the space of a/j3 vs a, as illustrated in Fig. 3 for the fit with (29). The 
valley is described by the equation a ~ (21.5 + 3.1 a/ (3) MeV. An analogous valley is found 
for the fit with A = 0, but the valley is then shifted by about 10% to higher values in a. The 
strong correlation may be understood in terms of the need to properly describe energies for 
the heaviest nuclei with largest asymmetries requiring a/(l + A^ 1 / 3 a/ (3) ~ a a , for A ~ 200. 

Along the valley, the average deviation minimizes at a/ (3 ~ 1.7, with the value of 
\E — E exp \ = 1.97 MeV. However, the minimum along the valley is rather flat. To un- 
derstand the nature of the minimum, we carry out the following procedure. We carry out 
an optimization of all parameters in the energy formula at a fixed a/ (3. Then, keeping all 
parameters fixed including at/ (3 but excluding a, we invert the energy formula to deduce 
a from the mass of an individual nucleus. With this, we test locally the global fit. If the 
formula provides a proper description, then, on the average, we should get the same a from 
an individual mass as from the global fit. Notably, for the inversion to be significant, the 
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asymmetry should be large and we demand that \N — Z\/A > 0.2. Sample results are 
displayed in Fig. 4, as a function of A. 

First to consider are the results obtained at a/ (3 — 0, i.e. ignoring the surface symmetry 
energy, as in the textbook mass formulas. Represented by the diamonds in Fig. 4, the 
results for a (i.e. a a in this case) group coarsely into a line. At large A, the line from 
symbols oscillates, exhibiting the shell structure, around the value of a from the global 
fit, represented by a straight line. However, at A < 50 the values from a local inversion 
systematically drop below a a from the global fit. Such a drop for a a at low A is expected 
from a competition between the volume and surface symmetry energies (28). If we raise the 
parameter ratio, such as to the value of a/ (3 — 2, we can get the locally restored values of a 
to oscillate around the global fit value down to very low A. Regarding the restored values 
of a, we may note in Fig. 4 that they locally exhibit less scatter for the fit with (29) than 
with A = 0, suggesting that the first fit grasps the physics better. Finally, if we assume a 
too high value of the parameter ratio, such as a/ (3 = 6, we get the restored a to rise up at 
low A in an opposite tendency to that for a/ (3 — 0. 

The gain in the description of low- A nuclei when using a finite a/f3 can be further 
appreciated by looking directly at magnitudes of the average fit deviations in the low- A 
region. Thus, if we employ a/ (3 = in (28) in the global fit with (29), then for nuclei in the 
region of 5 < A < 50 and \N-Z\/A > 0.2, we get \E - E ex P\ = 7.4 MeV, far in excess of the 
average deviation of 2.37 MeV for all nuclei. On the other hand, for a/ (3 ~ 1.7, we get an 
average deviation of \E - E ex P\ = 2.3 MeV in the region of 5 < A < 50 and \N-Z\/A > 0.2, 
close to the average deviation of 1.97 MeV for all nuclei. 

While, at the lower A, the use of an intermediate value of a/ (3 improves the performance 
of the energy formula, at the highest A the separate sensitivity to the ratio and a is quite 
reduced. With most of the measured nuclei of high anisotropies at the high A, the reduction 
in the deviation between the formula and the experimental energies, when averaged over all 
A, is just ~ 0.40 MeV per nucleus for the absolute minimum in Fig. 3, compared to the 
optimal case of a/f3 — 0. Considering in particular Fig. 4, we need to decide how accurate 
our phenomenological theory should be in describing the nuclei. While the superiority of the 
a//3 — 2 case over the other cases in Fig. 4 is quite apparent, cases that might be less than a 
third of the way in-between a/ (3 — 2 and the two other would not be clearly superior, given 
the scatter of the points and the lack of shell-effect description in the theory. With this, 
we consider the locus of points in the (a, a/ (3) plane corresponding the deviation reduced 
by a factor of 1/3 compared to the optimal a/ (3 — point (i.e. 0.13 MeV in Fig. 3), as the 
approximate boundary of the parameter range favored by the energy formula. Have we had 
a theory that could e.g. predict the shell effects, we might have been able to impose a more 
stringent limit on the parameters. For reference, we further show in Fig. 3 the boundary for 
the 0.40 MeV deviation, passing through the optimal a/(3 = point and in the vicinity of 
the optimal a/ (3 — 6 point. 

On one hand, the number of parameters here, in accessing the symmetry energy from the 
binding-energy data, represents a bare minimum, with the results necessarily being some- 
what coarse. On the other hand, there can be advantages to the fits with a low parameter- 
number, besides the obvious immediate accessibility of the analysis. While coarse, the low 
parameter-number fits tend to be stable and faithful, without the optimization being thrown 
off by obscure features of data or by auxiliary theoretical inputs. As an example of dangers 
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encountered when increasing the parameter number, when making the surface diffuseness d 
in (29) a variable parameter [13], we find that the average deviation from measured energies 
minimizes for sensible parameter values. However, that minimum is found to be very flat 
when the diffuseness and symmetry parameters are varied simultaneously, extending far out 
into the region of large surface-to-volume parameter ratios and large diffuseness. Such be- 
havior of a minimum marginalizes the significance of optimal values for a fit, given that, even 
when experimental errors are small, the theoretical uncertainties need to be incorporated 
[34]. The situation can be remedied by bringing in external knowledge [35], such as from 
electron scattering in the case in question. Such issues may be obvious for experimentalists, 
but are not routinely appreciated by theorists. 



Interest in asymmetry skins has surged in the context of the general interest in asymmetry 
effects, stemming from the availability of radioactive beams. Expectations were raised that 
an analysis of the asymmetry skins could permit to extrapolate the nuclear matter properties 
to the neutron matter [4-7] . Different features of effective interactions have been tested with 
regard to their bearing on the thickness of the skins. In several works [36] , a direct connection 
between the neutron skin in the spherical 208 Pb nucleus and the neutron star structure was 
made. Primarily, in recent times, the skins have been theoretically investigated within the 
microscopic structure calculations [37-39]. 

Unfortunately, experimental extraction of the neutron distributions has been difficult. 
Different methods have been employed, such as analyses of nucleon scattering [20,21,26,27] 
(including the analysis of polarized proton scattering and of neutron and proton scattering) 
or analyses combining information from strongly interacting (p, 7r ± , a, 12 C) and electromag- 
netic probes (e scattering, isotopic transition shifts) [22-25] . Recently, it has been proposed 
[40] to utilize parity violations in electron scattering to determine neutron rms radii, con- 
tributing to the surge of interest in the asymmetry skins. 



To assess the thickness of the asymmetry skin for a given surface excess (25), let us 
consider the sharp sphere radii for neutrons R n and protons R p . These radii correspond to 
volumes with interior densities such as for the bulk matter, i.e. 



In the last step, we use the fact that the intrinsic consistency of our considerations requires 
that we can tell the surface from the interior and, thus, must have \R n — R p \ C i? or, 
equivalently, \N S — Z s \ <C A We next get from (31) 



III. ASYMMETRY SKINS 



A. Interest in Skins 



B. Skin Size 




(31) 



Rn-R = Ns_ 

R 3N 



(32) 
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With an analogous result for R p , we find for the difference of the sharp-sphere radii: 

R n -R P = A (N 3 - Z s ) J_ N-Z-^ZAW 
R 6NZ 6NZ 1 + &A 1 / 3 



(33) 



where we inserted the result from (25). We see that the difference between the neutron and 
proton radii is primarily linear in the asymmetry and measures the symmetry coefficient 
ratio a/ 13. For a surface with no capacitance for the asymmetry, characterized by a/ (3 = 0, 
the asymmetry skin should vanish. If one were trying to compare the result (33) with the 
results in [14], one should note that there, in derivations, the approximation of iV ps Z ps A/2 
is utilized. 

Regarding (33), the results of measurements and of microscopic calculations for the skins 
are usually, however, not expressed in terms of the neutron-proton difference of sharp-sphere 
radii but rather the difference of rms radii. The difference between the rms compared to the 
sharp-sphere radii, aside from the standard reduction factor (dropping out from the ratio 
of radii) and from a correction for diffuseness, brings in a Coulomb correction additional to 
that already in (33). That additional correction stems from a polarization of the nuclear 
interior by the Coulomb forces. The polarization is schematically indicated in Fig. 2. In 
order to assess the polarization effect, within the net nuclear energy we represent the interior 
symmetry and Coulomb energies in integral forms. We expand these energies, respectively, 
to the second and first order in the deviation of asymmetry from uniformity: 

E = E' + yJdr (5p a ) 2 ~\J dr $(r) 5p a . (34) 

Here, $ is the Coulomb potential, p a is the neutron-proton density difference, p a = p n — 
p p , and we exploit the relation for the density changes Sp p = —5p a /2 following from the 
requirement p n +Pp = Po- The energy E' is the remaining portion of the net energy minimized 
with respect to the global partition of asymmetry between the interior and surface. The 
minimization of (34) with respect to 5p a , under the constraint of f dr5p a = 0, yields 



8a A 1 / 3 

where the overline indicates an average over the interior volume and where we utilize the 
sharp-sphere potential $. On calculating contribution of the deviation of asymmetry (35) 
from uniformity, to the difference of rms neutron and proton radii, we find 

(r 2 )V2 168a N ' [ } 

Notably, the contribution from the polarization (35) to the net energy is practically ignorable. 

On combining (33) with (36) and on incorporating a correction for the surface diffuseness 
into (30), we get for the difference of the neutron and proton rms radii: 

(r 2 )^ - {r 2 )T A N-Z a c A 5/3 f + g d (d n - d p ) 

{r 2 ) 1/2 6NZ1+PA 1 / 3 168a N 1 + £ A 1 / 3 (r 2 ) ' 1 ' 
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Calculations [4] appear to indicate that, for lighter nuclei, the difference in the neutron and 
proton diffuseness d n — d p is of the second order in the difference of sharp-sphere radii (33). 
For heavy nuclei, there may be a Coulomb contribution to the difference in the diffuseness. 
Overall, in any case, the correction for diffuseness is of a lower order in A than the leading 
term in the skin thickness. 

In Fig. 5, we compare the results of skin measurements in Na isotopes [25] to the results 
from (37) for different assumed values of a//3, with the difference in diffuseness set to zero 
and with a made to follow the binding-energy correlation valley in Fig. 3. The data rather 
clearly rule out a/ (3 < 1, but it is difficult, so far, to determine the energy parameter ratio 
with a much better precision. Further comparisons to data are shown in Figs. 6 and 7, for 
Sn and Pb isotopes, respectively. Overall, the displayed data appear to favor the ratios in 
the vicinity of a//3 ~ 3. As the charge number increases, one can see in Figs. 5-7, from 
the a/f3 = results, the increasing importance of Coulomb corrections to the skin size. 
For realistic ratios a/ (3 < 3, we find that the Coulomb forces contribute about 0.07 fm of 
reduction to the skin of 208 Pb. This can be compared to the ~ 0.1 fm change in the skin 
size found by Furnstahl [7] for 208 Pb when switching off the Coulomb interactions within 
mean-field models. The interior polarization and pushing out of the proton surface by the 
Coulomb interactions contribute about evenly to the skin reduction in our estimates. 

It needs to be mentioned that, when different methods are employed in the skin extraction 
from measurements, such as in the case of Ca isotopes [20-23,27], the outcomes can easily 
differ by more than the claimed extraction errors, underscoring the difficulty in extractions. 
Results from a global fit of (37) to the available skin sizes obtained from data [20-27] are 
next shown in terms of contour lines (solid, slanted downward) in Fig. 3. As has been already 
discussed, the dependence of skin sizes on asymmetry primarily constrains the values of the 
ratio a/ (3. We provide the contours for two combinations of a standard statistical deviation 
and an assumed average theoretical error S of (37), i.e. x 2 — Xmin + 1 + J2k(^/ cr k) 2 where k 
is data index, and 5 = 0.02 or 0.04 fm. 



C. Skin Size vs Difference of Separation Energies 

A strong dependence of the sizes of asymmetry skins on the proton-neutron difference 
of separation energies was predicted within the mean-field calculations of Ref. [41]. The 
dependence has been observed experimentally [25] and it was referred to in other mean-field 
calculations [42,5]. As we shall next demonstrate, this type of dependence naturally emerges 
in our results. The dependence can be used to constrain directly the surface parameter f3. 

The proton-neutron separation energy difference, up to a factor and 0(v4 -1 ) effects, is 
nothing else but the asymmetry chemical potential fj, a . From the energy contributions (17) 
and (24), we get right away 

_ dE _ 2(3{N S -Z S ) _ 5 Z 
** a ~ d(N s - Z S ) ~ A^ 6 a ° AVs • [6 *> 

The difference of sharp-sphere radii (33) can be then written in terms of the chemical po- 
tential as 

Rn-R P = ^ (N z) _ Pa ^ 5/3 ■ 5 ac A*/* 

R 6NZ { S S) 12(3 NZ 72 (3 N ' K ' 
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Turning to the difference of the r.m.s. radii, with the polarization contribution (36), the 
diffuseness correction and with \i a = (S p — S n )/2, where S p and S n are the proton and 
neutron separation energies, we find 

(r 2 )j /2 - (r 2 )^ 2 _ S p - S n A 5/3 a c A 5/3 f 35 a \ 2 d(d n -d p ) 

(r 2 >V2 24/3 NZ 168a N \ 3 (3 A 1 / 3 J (r 2 ) ' 1 ' 

If not for the Coulomb correction, the surface asymmetry in (38) would have been pro- 
portional to the difference of separation energies, with the coefficient of proportionality 
inversely proportional to (3. The linearity with respect to the separation-energy difference 
is largely retained for the skin size in (40). (Principally, one might strive in Eq. (40) to 
eliminate the separate dependence on iV and Z, in favor of exclusively S n — S p and A, but 
the separate dependence on iV and Z appears only in correction terms, either associated 
with the Coulomb forces or geometry, so it is not crucial.) Into the results previously com- 
pared to data, i.e. the dependence of energy and of skin size on mass and charge numbers, 
the surface and volume symmetry coefficients entered in a combination. By contrast, the 
dependence of skin size on the separation energy difference, that we now examine, tests 
primarily the surface symmetry coefficient. The difference between the cases is easily under- 
stood following the capacitor analogy. Thus, the volume and surface coefficients characterize 
the surface and volume capacitances, respectively. Either the net energy or the ratio of the 
surface-to-net asymmetry depended on the net capacitance. However, in the capacitor anal- 
ogy, the separation-energy difference represents a voltage and the relation between surface 
charge and the voltage depends only on the surface capacitance (as a matter of capacitance 
definition). 

We now confront the result (40) with the data on skin sizes as a function of separation- 
energy difference for Na, Sn and Pb nuclei in Figs. 8, 9 and 10, respectively. The data 
appear to favor the values /? > 10 MeV. In fitting all data on skin sizes [20-27], now with 
Eq. (40), we again generate contour lines (solid, slanted upward) in the plane of a-a/f3 in 
Fig. 3, for x 2 = Xmin + 1 + J2k(^/ cr k) 2 where k is data index, and the systematic error equal 
either to 5 = 0.02 or 0.04 fm. Since the dependence primarily constraints f3, the constraint 
lines in the plane a-a/[3 are approximately straight lines passing through the origin of the 
coordinate set. 



D. Symmetry Parameter Values 

With the three dependencies fitted to the data, we now have three regions in the plane 
of Fig. 3, favored by the data, that intersect. Since the separation energies stem from the 
binding energies, one would expect that the intersection of the two regions derived from 
skin sizes would be well centered over the region favored by binding energies. Some offset 
may be from a bias due to the fact that the skin data primarily stem from nuclei near 
closed shells or due to the the fact that the surface diffuseness is ignored in the Coulomb 
corrections for skin size but is not ignored in the binding-energy formula. In any case, the 
offset suggests that we may need to use bias errors as large as 5 = 0.03 fm for Eqs. (37) and 
(40). Intersection of the three regions favored by data then yields the following limits on 
the parameter values: 2.0 < a/ j3 < 2.8 for the ratio, 27 MeV < a < 31 MeV for the volume 
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and 11 MeV < (3 < 14MeV for the surface. Within those limits, though, the parameters are 
correlated, cf. Fig. 3. 

It should be mentioned that values of the volume symmetry parameters a in some ad- 
vanced mass descriptions fall within our region, such as in [43] with J = 28 MeV and in 
[13] with their — k vo i b vo i = 29.3 MeV. However, in [43] the particular parameter value is set 
and in neither of those sources the uncertainties are actually assessed. Elsewhere, values 
outside of our range are followed [12] and in the mean-field structure calculations, based on 
Skyrme interactions or Lagrangians, a wide range of the symmetry volume parameters is 
employed [7], from 24.5 to 43.5 MeV. Wide parameter ranges are also employed in various 
phenomenological considerations such as for the nuclear liquid-gas phase-transition [44]. On 
the other hand, nuclear-matter calculations relying on realistic nucleon-nucleon interactions 
appear to yield volume symmetry energy values within our range, of 27-29 MeV depending 
on the interaction in [45] and 30.10 MeV in [46]. Those microscopic energies, however, are 
obtained as the difference between the neutron and symmetric nuclear-matter energies but 
this, normally, is an excellent approximation to the symmetry energy from N — Z expansion. 

Regarding the obtained values of the ratio of a//3, one may ask about a justification for 
the large- A expansion of the symmetry energy in mass formulas [12,13]. With the expansion 
amounting to an approximation in (28) of 1/(1 + x) = 1 — x + x 2 — . . . ~ 1 — x, where 
x — a/(P A 1 / 3 ), we may adopt a criterion that the expansion is justified when the correction 
to the effect of surface on symmetry energy, from ignoring higher-order terms, is less than 
20%. This requires x < 0.2 or A > (5a / (3) 3 > 1000 which is never satisfied. (Even 
when the next term in the expansion is included [13], meeting the criterion necessitates 
A > 5 3 / 2 (a//3) 3 = (90-250), depending on the a/(3 ratio.) In the droplet model fit [12] 
to binding energies and fission barriers, with nonexpanded symmetry energy, the value of 
the ratio a//3 is within our region, as represented by | ^ = 2.5, and so is the value of 
the surface parameter /3, represented by | Q = 13.0 MeV. However, the uncertainties are 
not produced there and, moreover, with an inclusion of the volume symmetry parameter, 
J = 32.7 MeV, the results are well outside of our correlation region. In earlier mass and 
barrier fits within the droplet model, the ratio itself and the surface parameter were much 
different, with | ^ = 4.9 and | Q = 7.5 MeV in [16], respectively. It would be interesting 
to find out the value of the surface symmetry coefficient for the advanced mass calculations 
in [43]. In [13], the volume-to-surface symmetry parameter ratios a/f3 for mass fits are 
—bsurf K sur f/(b vo i k vo i) = 1.31 and 1.35, but these results are for an expanded form of the 
symmetry energy, at two different levels of the expansion. 

In the next section, we employ our elementary results in the description of asymmetry 
oscillations in nuclei. 

IV. APPLICATION TO ASYMMETRY OSCILLATIONS 

A. Goldhaber- Teller and Steinwedel-Jensen Models of the Giant Dipole Resonance 

So far, we discussed elementary static properties of ground-state nuclei. Nuclear exci- 
tations may involve displacements of protons relative to neutrons, compared to the ground 
state, with prominent giant resonance structures in cross sections allowing for an interpreta- 
tion in terms of a collective motion involving a number of nucleons. The lowest in energy (or 
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frequency) is the giant dipole resonance (GDR) prominent in the low-energy photoabsorption 
[47] and inelastic electron-scattering [48] cross-sections. 

Two early models have been advanced for the GDR, under the names of Goldhaber and 
Teller (GT) [29] and Steinwedel and Jensen (SJ) [29,30]. In the GT model, the neutron 
and proton distributions are assumed to oscillate against each other as rigid entities, see 
Fig. 11. Since the potential energy associated with the displacement is proportional to the 
nuclear surface area, i.e. oc A 2 / 3 , while the moving mass is proportional to A, the resonance 
excitation energy in this model is proportional to A" 1 / 6 : Eqdr — ^ oc a/ A 2 / 3 /A = A~ l l G . 
In the SJ model, a standing wave of neutron vs proton displacement develops within the 
nuclear volume, satisfying the surface condition of vanishing asymmetry-flux, analogous to 
the closed-pipe condition in the standard example of oscillations for a continuous medium. 
As the wavelength is then in a definite constant ratio to the nuclear radius, the resonance 
energy in the model is proportional to A^ 1 ^ 3 : Ecdr = foc a /\ oc A" 1 / 3 , where c a is the speed 
of propagation for asymmetry perturbations in normal symmetric matter. 

In detail, within the SJ model [30,49], it is assumed that there is no net movement 
of matter, i.e. the net matter velocity v vanishes and the net nucleon density stays intact, 
consistent with the continuity equation: dp/dt = — V(p v). The relative density p a = p n — p p 
changes, but the state of matter is locally the same as within a ground-state nucleus, just 
at the new asymmetry. The vanishing net velocity implies a relation between the velocities 
for individual particles, v n and v p , and densities, following from the consideration of fluxes 
that need to meet the condition 

= pv = p„v n + p p v p . (41) 

This then implies that the flux in the asymmetry continuity equation may be expressed in 
terms of the relative velocity v„ p = v n — v p , 



dpa 
dt 



- V(p„ v n - p p v p ) = -V I - v np 1 ~ Vv np 



(42) 



where the approximate equality is valid in the lowest order in the disturbance. An Euler 
type equation for the relative velocity can be next used to close the equation above, 

d 2 „ 2 dp a „ 2c 2 

JT + V « P = * V//a = 7T^Vp a = a - Vp a , 43 

ot m* m* op a Po 



where m* is the effective mass and 



2 Po dp a 



° a m* dp a ■ l44j 

On differentiating Eq. (42) with respect to time and on substituting (43), we get the wave 
equation for the asymmetry density 

d 2 4NZ „ 2 

dt 2pa = ^ ' a Pa - ( } 

Under spherical symmetry, the normal vibrations are looked for, in the standard way, 
in terms of the deviation from ground-state density expressed as the product of a spherical 
Bessel function jg and of a spherical harmonic Y^ m , 
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5 p a = D v sin (ut) je(qr) Y £m (Q) , 



(46) 



where, following (45), the wavenumber q and the frequency u are related by 

u = — - — c a q . (47) 

The radial velocity component, from (43) and (46), is 

2D V A 



i ■;, , = = c„ cos (ut)j' e (qr) V, „,(<>) . (48) 



The condition that the normal component of flux vanishes at the surface [28], v^ p (R, Q) = 0, 
yields the wavenumber values 

qen = § , (49) 

where n orders the values and j' £n represents the n's zero of the spherical Bessel function 
derivative, = 0. For the resonance energy, we get, given (47), 

t, ln — Ruin — — — — nc a qi n — — , (bU) 

where, with (44) and (16), 




(51) 

The GDR is characterized by £ = 1 and n = I. With j' u ~ 2.08, the GDR energy in the SJ 
model is then, finally, 

2(NZ) 1 / 2 f2a\ l/2 2.08 k ( 2a \ 1/2 2.08^ . , 

E GDR= \ A L — -— « — — TT7T- (52) 



A 4 / 3 \m* J tq \m* J tq A 1 / 3 

Within the GT model of the GDR, the changes in asymmetry occur only within the 
nuclear surface and in the SJ model only within the interior. Experimental dependence 
of GDR energy on mass number, displayed in Fig. 12, is intermediate [47] between that 
characteristic for the GT and SJ models, suggesting that both the volume and surface 
participate in the oscillations. To explain the mass dependence of the GDR energy, in 
Ref. [50] the oscillation in the SJ model was treated as an elementary mode of nuclear 
oscillation coupled to a second GT mode of oscillation, within the same nucleus. As the 
resonance frequencies for the two models, GT and SJ, shifted with a changing nuclear mass, 
the lower-frequency normal mode for the combination of oscillators changed from a mode 
with a predominantly GT to SJ content. However, this casts the nuclear system into one 
with only two degrees of freedom. The physical situation is rather that of a continuous 
system for which, though, the surface can generally accept the asymmetry flux, making the 
original SJ boundary condition excessively restrictive. 
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B. Description of Asymmetry Oscillations 



We now explore the modification and consequences of the surface boundary condition 
for the asymmetry oscillations. Given that the surface can store an asymmetry excess, we 
replace the SJ condition of a vanishing flux with the condition for the change in surface 
asymmetry density n a , 

d r r 2NZp r 

^n a = p n v n - p p v p = A2 v np , (53) 

cf. (42). As the chemical potential for the surface needs to be the same as within the 
adjacent interior, this can be further rewritten as a condition involving the interior chemical 
potential, or involving the volume density next to the surface, 

2NZp r d_ dUa_ Opa = dt}a_ Opa dpa , .s 

A 2 ~ dt^' dp a dt ~ dp a d Pa dt ' 1 J 

When the surface does not accept the excess, i.e. dn a /dp a = 0, we get the regular SJ 
condition of a vanishing flux. 

Looking for normal modes with (46) and 

8n a = D s sm(ut) Y em (Q) , (55) 

we find from the requirement of equality for chemical potentials at the surface, in accordance 
with local surface and volume densities of symmetry energy, that 

Ds = ^ JU - 3 D v Rj e (qR). (56) 

Moreover, with (48) and (47), we can obtain from either (53) or (54) the condition for the 
wavevector of oscillation: 

3/541/3 

qRj e (qR) = ^ j' e (qR) . (57) 

a 

As the ratio a/ (3 A 1 / 3 grows from small to large values, the relative magnitude of the 
surface amplitude D s in (56) grows. Moreover, it is apparent that the solutions to the 
wavenumber equation (57) generally shift then from j' in / 'R to je n /R, where ji n represents 
the n'th zero of the spherical Bessel function, je{je n ) = 0. In terms of the standard example 
of oscillations for a continuous medium, the boundary condition shifts from one representing 
a closed pipe to one representing an open pipe. Coping with the n — 1, £ > solution to 
the wavenumber equation at q — > 0, as (3 A 1 / 3 /a — > 0, however, requires a careful expansion. 

We concentrate next on the £ = 1 GDR case. For a/ f3 A 1 ' 3 -> 0, the surface amplitude 
in (56) tends to zero and we arrive at the SJ limit with the GDR energy given by (52). For 
(3 A 1 / 3 /a — > 0, we need to consider the Bessel function in (57) at small values of its argument 
where ji(x) ~ (x/3)(l — x 2 /10). We then find for the wavenumber 

^© V2 ^ (58) 
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and for the resonance energy 



2VNZ _ 2(NZy/ 2 {Qf3\ 1/2 h /6/A 1/2 h 

EGD " = ^ = —^ h ^ = —^\m-) V^W) ^W*' (59) 

In this limit, the resonance energy becomes dependent only on the surface symmetry param- 
eter (5 and not on the volume parameter a. Moreover, the volume amplitude gets suppressed 
relative to the surface amplitude in (56). This is clearly the GT limit of vibration. 

Figure 12 shows next fits to the data [47] on mass dependence of the GDR energy, with 
c a treated as a fit parameter and with q u obtained from (57) for different values of a//3. 
We see that a realistic ratio of a//3 — 2.5 yields a much better description of the data than 
the SJ ratio of a//3 — 0. However, a still larger ratio than 2.5 would provide an even better 
description. It is possible that other factors than the finite capacitance of the surface for 
asymmetry play a role here, such as the local deviations of a system from the ground state 
[51]. For the best fit with a/f3 = 2.5, we get c a = 0.28c, corresponding to m — 0.84m*. 

One may ask for what values of A and realistic a/f3 the vibrations turn from more GT 
to more S J like. As a criterion, one may use the fraction of potential energy of the vibration 
attributable to the surface. With this criterion, one finds that the transition takes place 
at pAH 3 /a = 3/2, and, thus, for the range of 2.0 < a/0 < 2.8, at 30 < A tr < 75. The 
relatively high transition masses can be attributed to the node of the vibration at the center 
of a nucleus. 



C. Transition Densities 

The profile of a standing wave in an oscillation represents, at the quantal level, the tran- 
sition density for the resonance. This density is probed in inelastic electron scattering [48]. 
To our knowledge, no inversion of data for the GDR transition density has been ever made, 
but the data were tested against the GT and SJ type densities, with a conclusion that the 
GDR (and also £ > 1) densities are in-between the two limits [48]. Even in theoretical works, 
the transition densities have not been much quoted. 

We write the transition density as 



Pi(r) = pi(r) Y tm , 



(60) 



where pi(r) includes both the volume and surface contributions. Interpreting the volume 
results in terms of long-wavelength modulations of the particle densities and the surface 
results in terms of the modulations of radii, we have for the density [48] 



Pi(T) = — 
Po 



dp 

Dv Je{qr) p{r) - D s — 



(61) 



where p{r) is the ground-state density distribution. With (56), we can further rewrite this 

as 



Pi(r) = —jt(qr) 
Po 



p(r) 



a 



dp 



(62) 
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Establishing a contact with microscopic theory, for 40 Ca in Fig. 13, we compare our GDR 
transition density (62), with (30), to the transition density from the microscopic calculations 
of Ref . [52] , that include the effects of 2p-2h excitations and of ground-state correlations. The 
normalization in (62) is obtained from the Thomas-Reiche-Kuhn sum rule; for illustrative 
purposes, we take m* = m as in [52]. The calculations [52] employ a phenomenological 
Landau- Migdal type interaction in use in describing different giant resonances. We note 
that a best contact with our results is made for a/ (3 > 2.5. 



V. DENSITY DEPENDENCE OF SYMMETRY ENERGY 



Several works tried to identify links between the asymmetry skins and bulk asymmetry 
properties of nuclear matter, within microscopic theory [4-7] . Having explored a link between 
the skins and the symmetry parameters, we will now try to establish a connection between 
the parameters and the density dependence of symmetry energy, starting with the Thomas- 
Fermi theory. 



A. Symmetry Parameters in the Thomas-Fermi Theory 

We shall consider a nuclear energy functional expanded to the second order in asymmetry: 

E = E + J drpS(p)(^j =E + E a . (63) 

Here we assume that the expansion term in asymmetry is local, but this will not be essential 
for much of the discussion that follows. The function S is the density-dependent symmetry 
energy per nucleon and S(po) = ol. The energy E is generally nonlocal but independent 
of asymmetry as the zeroth-order term in the expansion; since the Coulomb effects have 
been explicitly considered in Sees. II and III, we will disregard them for the discussion that 
follows. Regarding the role of asymmetry in any nonlocality, there had been indications 
[4,7,53] that it has a negligible impact on the properties of nuclei; the apparent success in 
systematizing the asymmetry skins in terms of the nuclear-matter properties [4-7] lends, in 
itself, a support to the weak impact. 

A simple example of the functional (63) is such as employed in the initializations of 
transport calculations for central reactions [54]: 

E = [dr(le F p + ^ + 7 ^—) + ^[dr(Vp)\ (64) 
J V 5 Po ( (7 + 1 )Po/ 2 Po J 



S '+ C -l + ^, (65) 
3 Po Po 



with the Fermi energy ep = (3tt 2 /2) 2 / 3 h 2 p 2 / 3 /2m. The constants A, B and a are fixed by 
the requirement that the functional yields the energy minimum of E/A = ay ~ —16 MeV 
at po in symmetric nuclear matter, characterized by a specific incompressibility constant K. 
E.g. for K = 260 MeV, we get A = -180.9 MeV, B = 128.1 MeV and a = 1.446. The 
constant V provides a size scale for the surface diffuseness and an adequate surface profile 



22 



is obtained [54] for T> k> 20.4 MeVfm 2 . As to the constants in S, we have the relation 
a = €f(po)/3 + C + T and, otherwise, different combinations of the constants C and V yield 
different values of (5 as will be discussed. 

The minimization of the energy (63) under the conditions of fixed particle numbers 

A = j drp and N-Z = J dr p a (66) 

produces the Thomas-Fermi equations: 

5E 2 d S 

»=li + p «T p -p> (67) 

p a = 2p a ^-. (68) 

At first, those equations will just serve as a background illustration for the claims made in 
the discussion. 

On the basis of (63), we would like to find out changes in the nuclear ground-state energy, 
as the particle asymmetry changes. First, we consider a symmetric system and find p(r) 
that minimizes the energy E = E . With the energy E being quadratic in asymmetry, 
the net density minimizing the energy functional at a finite asymmetry will differ from the 
density for the symmetric system only to within a second order in asymmetry. With this, to 
within the second order in asymmetry, the ground-state E a in (63) may be calculated with 
p(r) obtained for symmetric matter. Moreover, since the symmetric-system p(r) minimizes 
E , any changes to it, of the second order in asymmetry, will produce only quartic terms of 
change in E . In conclusion, if we are interested in the change in the ground-state energy 
to the second order in asymmetry, we just need to consider the E a term in (63) and we can 
use there the density p(r) from the energy Eq minimized for a symmetric system. In terms 
of the Thomas-Fermi equations, this means that we can uncouple (67) by putting p a = 
there. 

From the second of the Thomas-Fermi equations, we have in the ground-state system 



and from (63) 



* = T§ (69) 



B« = f (70) 



Within the analogy with a capacitor, the integral j dr(p/2S) represents the system ca- 
pacitance, independent of the deposited charge. The volume capacitance is A/ 2a and the 
reminder is the surface capacitance: 

A2/3 s 1 r, ( i i \ 

- drp --- . (71) 



2(3 2nrlf3 2 J ' \S a 

As we will not be pursuing here curvature effects, we just turn to the surface capacitance for 
a semiinfinite matter. With x a coordinate along the direction perpendicular to the surface 
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and on factoring out the density dependence from the symmetry energy with S — a s(p/po), 
we find from (71) 



Thus, for (63), given that the density profile for the surface is constrained, the ratio of 
the symmetry parameters a/ (5 is a measure for the density dependence of the symmetry 
energy S. We note that a density-independent symmetry energy, characterized by s = 1, 
yields a/ (3 = 0. The falling density at the surface produces in this case no extra capacitance 
per nucleon compared to the interior. The faster the drop of s with density, the larger the 
extra capacitance, exhibited in a larger a//3 ratio. 

We next test the validity of the considerations from Sees. II and III within a Thomas 
Fermi model based on Eqs. (64) and (65), with a minor alteration, however, that the different 
neutron and proton Fermi energies are treated explicitly for asymmetric nuclei. On one hand, 
we solve the differential equation [54] from (67) for the semiinfinite symmetric matter and 
we vary C and T to get various combinations of a and a//3 from (72), requiring that they 
follow the binding-energy correlation from Fig. 3. Given the a and (5 values, we make a 
prediction for the size of the asymmetry skin of a specific nucleus, following Eq. (37) with 
d n — d p = 0. On the other hand, we solve the Thomas- Fermi equations directly for a nucleus, 
with Coulomb forces included, and calculate the difference of rms radii. A comparison of the 
results for Na isotopes is shown in Fig. 14. The results are both compared when changing 
the isotope mass, while keeping the symmetry parameters fixed, and for one isotope when 
changing the symmetry parameters. As is seen, the agreement of the estimate (37) with 
the actual skin size is rather remarkable. For heavy nuclei, however, a small systematic 
difference develops between the actual and estimated skin sizes, that largely vanishes when 
the Coulomb forces are switched off. The size of the discrepancy, ~ 0.03 fm for Pb, for which 
we have no understanding, is compatible with the possible difference between our estimate 
of the Coulomb effects in Eq. (37) and that of Furnstahl [7] (already brought up in Sec. III). 

A very recent, more involved, analysis of surface symmetry energy and skin sizes, within 
the Thomas- Fermi theory, has been carried out by Bodmer and Usmani in Ref. [53], with 
some conclusions comparable to those that emerge here. In detailed comparisons of their 
results to ours, it should be noted that their surface symmetry energy and parameter are 
differently defined than ours for a finite nucleus and e.g. include Coulomb contributions. 



To gain a perspective on the recent microscopic results on skin systematics in the litera- 
ture [5-7] , we further confront our results from the Thomas- Fermi model, to those obtained 
in the variety of mean-field models explored there. In Fig. 15, the skins from the mean-field 
calculations for 208 Pb [7] (symbols) are presented vs the reduction factor s for the symmetry 
energy per nucleon at po/2 (left panel), vs the scaled derivative of the energy at normal 
density po = p^dS/dpo\ po (center panel) and vs the volume symmetry parameter a (right 
panel). The skin dependence on p and on a was explored in [7]. A similar derivative to 
Po was utilized in [5,6]. Moreover in Fig. 15, the results from the Thomas- Fermi model 




(72) 



B. Constraints on Density Dependence of Symmetry Energy 
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are presented (lines) for the symmetry parameters changing along the correlation valley in 
Fig. 4. It is apparent in each panel of Fig. 15, that the Thomas-Fermi model produces 
general trends such as found in the mean-field calculations. However, in the three panels, 
the represented dependencies for skins in the Thomas-Fermi model are of a quite different 
nature. Thus, when a just provides a normalization for the symmetry energy, it has a very 
limited bearing on the skins in (37), merely through a rescaling of the Coulomb corrections 
(see also [53]). However, when a is correlated with a//3 through the binding-energy require- 
ments, the changes in a force changes in a/j3 and, thus, produce significant changes in the 
skin in the rightmost panel of Fig. 15. In the center panel for the Thomas- Fermi model, we 
observe combined indirect and direct (albeit at the edge of the relevant density interval) ef- 
fects of the abscissa on the skin. Only in the leftmost panel for the Thomas- Fermi model, we 
observe the effect of a direct physical connection between the abscissa and the skin, through 
Eqs. (37) and (72). The three panels illustrate potential dangers of scouting for a physical 
relation [5-7] looking at numerical results within models that are necessarily constricted. 

For a given effective interaction, the mean-field models likely yield more realistic results 
in the surface region than the Thomas-Fermi model. However, there is no reason to believe 
that a strong primary sensitivity to the density dependence of s disappears, for the symmetry 
parameter ratio a/f3, in the mean- field models. Using the results [7] for the 208 Pb skins in 
those models, we can assess a/ (3 for the models, from (37), and we can produce constraints 
on the low-density s from the constraints on a//3. An alternative would be to use the data 
on 208 Pb skins directly to constrain s in Fig. 15. In Fig. 16, the values of s at po/2 are 
shown, for the mean- field models, vs a//3 from inverting (37). In the inversion, basing here 
the results on one heavy nucleus, we augment the Coulomb correction in (37) by 0.03 fm; 
an ad hoc procedure suitable for all nuclei might be to amplify the Coulomb correction by 
a factor of 1.4. Given the variation of the results in Fig. 16, for the acceptable range of 
2-0 < a./j3 < 2.8, we deduce the range of the symmetry-energy reduction factor at po/2 of 
0-57 < s < 0.83. Compared to the previous unbiased limit [7] on p^, of lMeV/fm 3 < p^ < 
4MeV/fm 3 , our limit on (1 — s) represents a nearly twofold improvement in constraining the 
low-density dependence of symmetry energy. The s-shift is for the fairness in comparison; 
note that the majority of theoretical calculations in Fig. 16 suggests an even narrower range 
for s than above, of 0.67 < s < 0.83. Concerning the microscopic nuclear-matter calculations 
with realistic nucleon-nucleon interactions, minimally extrapolating underneath the density 
range in [45], the reduction in symmetry energy at po/2, for the variety of interactions there, 
is within our narrow reduction range, at 0.67-0.73. However, in the calculations of Ref. [46], 
the reduction in symmetry energy at p /2 is only within our wider range, at 0.59. 

VI. CONCLUSIONS 

Binding in the symmetric nuclear matter can be realistically assessed by anyone with 
a minimal experience, using the textbook mass formula and a sample of nuclear masses. 
We hope to have demonstrated that the information on symmetry energy can be accessed 
with a nearly comparable ease, using a minimally modified formula. Elementary but general 
considerations, for the nucleus as a macroscopic system, indicate that the dependence of 
surface tension on nuclear asymmetry necessarily implies the emergence of an asymmetry 
skin. In the lowest order in asymmetry, the surface tension and surface energy need to be 
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quadratic in the asymmetry. With analogous behavior of the volume energy, the surface 
and volume symmetry energies combine as energies of two connected capacitors with the 
volume capacitance proportional to the volume and the surface capacitance proportional to 
the area. 

In minimizing the net energy, the net asymmetry partitions itself into the surface and 
volume asymmetry in proportion to the capacitances. As nuclear size increases, the surface 
portion of the asymmetry decreases, as surface decreases relative to the volume. The minimal 
net symmetry energy is given in terms of the net asymmetry squared divided by the sum of 
the volume and surface capacitances. 

In analyzing the binding energies, we demonstrate that the surface symmetry contribu- 
tion is required to describe the energies of light asymmetric nuclei at a level similar to that 
typical for other nuclei. Overall, however, the binding energies are better in constraining 
a combination of the coefficients for the surface and volume symmetry energies, than in 
constraining the individual coefficients, at least at the level of our simple formula. The 
specific coefficient combination represents the symmetry coefficient in the standard mass 
formula. On the other hand, the dependencies of skin sizes on asymmetry and on the 
neutron-proton separation-energy difference primarily constrain, respectively, the ratio of 
the symmetry volume to surface parameters and the surface symmetry parameter. Using 
the three types of constraints, we establish limits on the volume symmetry parameter of 
27MeV < a < 31 MeV, on the surface symmetry parameter of 11 MeV < (3 < 14MeV, and 
on the parameter ratio, of 2.0 < a/ (3 < 2.8. Within those limits, however, the parameters 
are correlated. In analyzing the skin sizes, we provide straightforward interpretations for 
regularities observed both in measurements and in microscopic structure calculations. 

For reference, when taking the parameter ratio in the middle of the permitted range, 
at a/j3 ~ 2.4, the sample best-fit parameter values in describing the binding energies 
with (28) and (29) are a v = 15.4808 MeV, a s = 18.4359 MeV, a c = 0.696424 MeV, 
a P = 9.82356 MeV, a = 29.1048 MeV and (3 = 12.1270 MeV. When ignoring the sur- 
face diffuseness in putting A = in (28), the best-fit value are ay = 15.6163 MeV, 
a s = 17.9878 MeV, a c = 0.692791 MeV, a P = 10.8714 MeV, a = 32.6655 MeV and 
P = 13.6106 MeV. 

We point out that the additional capacitance for asymmetry provided by the nuclear 
surface is associated with the drop of symmetry energy per nucleon with density. The ratio 
of the volume to surface symmetry coefficients can be used to assess the pace of that drop. 
We test the analytic result for the skin size within the Thomas-Fermi approximation and 
we use the results from microscopic mean-field calculations to constrain the drop in the 
symmetry energy at half of the normal density, s = S(p /2)/S(p ), finding 0.57 < s < 0.83. 

For the collective asymmetry oscillations, we derive a simple boundary condition, fol- 
lowing a consideration of the surface fed by the asymmetry flux from the nuclear interior. 
The nature of the condition depends on the nuclear mass, as the relative capacitance of the 
surface for asymmetry changes with the system size. As the size increases, the nature of the 
collective oscillation changes from a surface to a volume mode. 
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FIG. 1. Illustration for the construction of surface quantities following Gibbs [18]. The top 
shows the variation of the density / for some quantity F within the system. The bottom shows 
an idealized reference system, where the density changes in a discontinuous manner, right at the 
surface, between the density values fy characteristic for the system away from the interface. The 
surface quantities Fg are defined as the difference between F, for a macroscopic volume containing 
the interface, and F^, for the same volume in the reference system, Fs = F — F^. 
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R r 

FIG. 2. In a two-component system, the surface-attributable particle-asymmetry can arise 
even when the surface is constructed to have no net particle number. For a spherical ground-state 
nucleus, this corresponds to the situation when the sharp-sphere equivalent proton R p and neutron 
R n radii differ from each other and bracket the matter sharp-sphere radius R. The dashed lines 
in the top part of the figure illustrate the polarization of a nuclear interior, in the form of a 
sharp-edged sphere of radius R, by the Coulomb interactions. 
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a (MeV) 

FIG. 3. Results from optimizations of the volume parameter a and the volume-to-surface pa- 
rameter ratio a/ (5. The thick solid diagonal line shows the optimal values of a at fixed a/ (3 when 
fitting the nuclear masses. The oval contour lines show the loci of parameter combinations that 
yield average absolute deviations from the measured nuclear masses larger by 0.13 and 0.40 MeV, 
respectively, than a best-fit parameter combination. The rectangular contour lines indicate the loci 
of parameter combinations at constant x 2 when fitting the asymmetry skins deduced from data, 
either following Eq. (37) or (40). Only the solid more horizontal portions of the skin-thickness 
contours are meaningful - the vertical dashed portions serve only to guide the eye. The loci are 
for the x 2 - varues assuming a combination of the standard statistical deviation and of the average 
systematic error of a skin-size formula of 0.02 and 0.04 fm, respectively, with the solid positively 
sloped contours representing Eq. (40) and the negatively sloped representing Eq. (37). 



32 



> 



30 



20 



i — i — i — i — I i i i i I i i i 1 1 1 1 1 1 [ i ii 1 1 1 mi |Mii|ini|iiin 1 — i — i — i — I i i i. 




a/0=O 























Mil 


-H-H- 


4+++ 


4+tt 


tttt 


tttt 


HH 


tttt 


m 


1 1 1 1 | 



+++ 




a e (A,Z) 




a//S=0 



J I I I I i i i i I i i i i I i 1 1 1 ■ ■ ■ ■ ■ ■ ■■ ■■■■■■■■■■■■■■■■■■ i i i i I i i i 



10 20 30 50 70 100 200 

A 

FIG. 4. Values of the volume symmetry parameter a from inverting the energy formula at 
different values of a/f3, for individual nuclei with \N — Z\/A > 0.2, as a function of A. The top 
panel shows results for the energy formula with A = and the bottom for the formula with A given 
by Eq. (29). The different symbols represent results from the inversions for a/f3 = 6 (squares), 2 
(circles) and (diamonds), respectively. The horizontal lines represent values of a from the global 
fits at the different a/ ft. 
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FIG. 5. Difference between the neutron and proton rms radii for Na isotopes as a function of 
the mass number. The symbols represent results from the data analysis Ref. [25]. The different 
lines represent results from Eq. (37) for the indicated values of a//3, d n — d p = and a taken from 
the binding-energy correlation valley in Fig. 3. 
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FIG. 6. Difference between the neutron and proton rms radii for Sn isotopes as a function of 
the mass number. The symbols represent results from the data analysis of Ref. [20]. The different 
lines represent results from Eq. (37) for the indicated values of a//3, d n — d p = and a taken from 
the binding-energy correlation valley in Fig. 3. 
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FIG. 7. Difference between the neutron and proton rms radii for Pb isotopes as a function of 
the mass number. The symbols represent results from the data analysis of Refs. [24] (circles), [20] 
(a diamond), [27] (a square) and [26] (a cross). The different lines represent results from Eq. (37) 
for the indicated values of a/ (3, d n — d p = and a taken from the binding-energy correlation valley 
in Fig. 3. 
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FIG. 8. Difference between the neutron and proton rms radii vs the difference between the 
proton and neutron separation energies for Na isotopes. The symbols represent data of Ref. [25]. 
The lines represent predictions of Eq. (40) for different indicated values of (5 (in MeV) with the 
values of a made to follow the binding-energy correlation in Fig. 3. 
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FIG. 9. Difference between the neutron and proton rms radii vs the difference between the 
proton and neutron separation energies for Sn isotopes. The symbols represent data of Ref. [20]. 
The lines represent predictions of Eq. (40) for different indicated values of (5 (in MeV) with the 
values of a made to follow the binding-energy correlation in Fig. 3. 



38 



t — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — r 




I i i i i I i i i i I i i i i I i i i i i i i i ' i i i i 

-15 -10 -5 5 10 

S p -S n (MeV) 

FIG. 10. Difference between the neutron and proton rms radii vs the difference between the 
proton and neutron separation energies for Pb isotopes. The symbols represent results from the 
data analysis of Refs. [24] (circles), [20] (a diamond), [27] (a square) and [26] (a cross). The lines 
represent predictions of Eq. (40) for different indicated values of j3 (in MeV) with the values of a 
made to follow the binding-energy correlation in Fig. 3. 
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FIG. 11. In the GT model [29] of the GDR, neutron and proton distributions are assumed to 
oscillate, unchanged with respect to their centers, against each other. In the SJ model [29,30] of 
the GDR, particle asymmetry oscillates within the nuclear interior, with the flux of asymmetry 
vanishing at the nuclear surface. 
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FIG. 12. Energy of giant dipole resonance as a function of mass number A. Symbols represent 
photoabsorption data from the compilation [47]. For mass numbers A > 50, we use the energies 
from the Lorentzian function fits from [47] . When a splitting characteristic for a deformed nucleus 
is present, we plot the equivalent energy for a spherical nucleus. For mass numbers A < 50, we 
plot the peak values for nuclei for which a clear isolated peak is present in the photoabsorption 
cross section [47]. The lines represent energies calculated from Eqdr = (2\/ (NZ)/A) hc a qu with 
c a being a fit parameter and with qu obtained from (57) at the different indicated values of a/ (5. 
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FIG. 13. GDR transition density pi for 40 Ca, as a function of the distance r from nuclear 
center. The lines marked with the respective values of parameter ratio a/ j3 represent the densities 
obtained from Eq. (62). The line marked GSC represents density from the microscopic calculations 
of Ref. [52], including effects of the 2p-2h excitations and ground-state correlations. 
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FIG. 14. Size of the asymmetry skin estimated with Eq. (37) vs size obtained directly from the 
calculated nucleon distributions, within the Thomas-Fermi calculations of sodium isotopes. The 
filled circles represent the results for a//3 = 3, while changing the isotope mass. The open circles 
represent the results for 30 Na, when changing the values of a./f3 and correlating the changes in a 
to follow the valley in Fig. 3. 
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FIG. 15. Size of the asymmetry skin in models vs the reduction factor s for the symmetry 
energy at po/2 (left panel), vs the scaled derivative [4,7,53] of the energy po = p^{dS/dp) at po 
(center panel) and vs the volume symmetry parameter a (right panel) . Symbols indicate results for 
the variety of mean- field models explored by Furnstahl [7]. Lines indicate Thomas- Fermi results, 
from Eq. (37) using (72). 
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FIG. 16. Reduction factor s for the symmetry energy per nucleon at po/2 vs the symmetry 
parameter ratio a/ (5 deduced for the mean- field models from the calculated [7] 208 Pb skins and 
Eq. (37). 
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